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ABSTRACT 



Thanks to their enormous energy release which allows to detect them up to very 
high redshift, Gamma Rays Bursts (GRBs) have recently attracted a lot of interest to 
probe the Hubble diagram (HD) deep into the matter dominated era and hence comple- 
ment Type la Supernoave (SNela). However, lacking a local GRBs sample, calibrating 
the scaling relations proposed as an equivalent to the Phillips law to standardize GRBs 
is not an easy task because of the need to estimate the GRBs luminosity distance in 
a model independent way. We consider here three different calibration methods based 
on the use of a fiducial ACDM model, on cosmographic parameters and on the local 
regression on SNela. We find that the calibration coefficients and the intrinsic scatter 
do not significantly depend on the adopted calibration procedure. We then investigate 
the evolution of these parameters with the redshift finding no statistically motivated 
improvement in the likelihood so that the no evolution assumption is actually a well 
founded working hypothesis. Under this assumption, we then consider possible system- 
atics effects on the HDs introduced by the calibration method, the averaging procedure 
and the homogeneity of the sample arguing against any significant bias. We neverthe- 
less stress that a larger GRBs sample with smaller uncertainties is needed to definitely 
conclude that the different systematics considered here have indeed a negligible impact 
on the HDs thus strengthening the use of GRBs as cosmological tools. 

Key words: gamma -rays burst : general - distance scale - cosmological parameters 



INTRODUCTION 



The observational evidences accumulated in these less 
than fifteen years, from the anisotropy and polariza- 
tion spectra of the cosmi c microw a ve ba c kground radi- 
ation ( CMBR, [dc Bcmar dis et al. I (|200(]| ): iBrown et all 
(|2009l ): iKomatsu et al. I (|201ll ). the large scale structure 
traced by galaxy re d shift surveys JPodelson et al. 1 20021: 
Percival et aTl |2002| : ISzalav et al. I l2003l: iHawkins et al 



the many efforts and the plethora of models proposed along 
these years. Ironically, the problem here is not the lack of a 
well established theory, but the presence of too many viable 
candidates, ranging fro m the classical co smological constant 
(Ca rroll ct al. . 1992: Sahni fc Starobinskv , ,20001) to scalar 



20031). the matter po wer spectrum (|Tegmark et al. I l2006l : 
Percival et al. I |2007| ) with the imprints of the Bar; 



■ ry 

onic Acoustic Oscilla tions (BAO. [Eisenstein et al. I (|2005h 
iPercival et al. ] (|2010l)l and the Hubble dia gr am of Typela 



fields (jPeebles fc Ratra I I2OO3I : ICopeland et al. I I2OO6I) and 
highe r order gravity theories ([Capozziello fc Francaviglial 



Noiiri fc Odintsov I l200a 



Super n ovae (SNela, Kowalski et al. I (|2008l ) ; iHicken et al. I 
(|2009l ): iKessler et al. I (|2009l )). definitely support the cos- 
mological picture of a spatially flat universe with a subcrit- 
ical matter content (fijv/ ~ 0.3) and undergoing a phase of 
accelerated expansion. While the observational scenario is 
now firmly established, the search for the motivating the- 
ory is, on the contrary, still in its infancy notwithstanding 



20081: ISotiriou fc Faraonil I2OIC. . . . , 

De Felice fc Tsuiikawa Il2010l ). all of them being more or less 
able to fit the available data. 

As often in science, adding further data is the best 
strategy to put order in the confusing abundance of the- 
oretical models. In particular, pushing the observed Hubble 
diagram to higher redshift would allow to trace the universe 
background evolution up to the transition regime from the 
dark energy driven speed up to the decelerated matter dom- 
inated phase. Moreover, being the distance modulus /i(z) 
log -linearly related to the luminosity distance and depend- 
ing this latter on the dark energy equation of state through 
a double integration, one has to go to large z in order to 
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discriminate among the prediction of diflerent models when 
these predict similar ijl{z) curves at lower redshift. Unfortu- 
nately, SNela are not ideally suited to this task with their 
present day Hubble diagram going back to Zmax ~ 1-4 and 
not extending further than z ~ 2 even for excellen t space 
based projects such as SNAP (|Aldering et al. 112003 1. 

Thanks to their enormous almost instantaneous en- 
ergy release, Gamma Ray Bursts (GRBs) stand out as 
ideal candidates to go deeper in redshift , the farthest one 
yet being at z = 8.3 (|Salvaterra et al."1 12009| ). The wide 
range spanned by their peak energy makes them every- 
thing but standard candles, but the existence of many ob- 
servationally motivated correlations between redshift de- 
pendent quantities and rest frame prop erties |Amati et al. I 
20081: iFenimore fc Ramirez -Ruiz Il2000l: iNorris et al. IbOOg 



Ghirlanda et al. II2004I : iLiang fc Zhang l^2005^ offers the m- 
triguing possibility of turning GRBs into standardizeable 
candles just as SNela. The use of these scaling relations al- 
lows to infer the GRB distance modulus with an error mainly 
depending on the correlation intrinsic scatter. Combining 
the estimates from different correlations, Schaefer (2007) 
first derived the GRBs Hubble diagram (hereafter, HD) for 
69 objects, while Cardone et al. (2009) used a different cali- 
bration method and add a further correlation to update the 
GRBs HD. Many attempts on using GRBs as cosmological 
tools have since then been performed (see, e.g., Firmani et 
al. 2006, Liang et al. 2009, Qi & Lu 2009, Izzo et al. 2009 and 
refs. therein) showing the potential of GRBs as cosmological 
probes. 

It is worth stressing that the possibility offered by GRBs 
to track the HD deep into the matter dominated epoch does 
not come for free. Two main problems are actually still to 
be fully addressed. First, lacking a local GRBs sample, all 
the above correlations have to be calibrated assuming a fidu- 
cial cosmological model to estimate the redshift dependent 
quantities. As a consequence, the so called circularity prob- 
lem comes out, that is to say one wants to use GRBs scaling 
relations to constrain the underlying cosmology, but needs 
the underlying cosmology to get the scaling relations. Differ- 
ent strategies have been proposed to break this circularity 
so that our first aim here is to investigate whether they 
are indeed viable solutions and to what extent the residual 
problem impact the HD derivation. 

A well behaved distance indicator should be not only 
visible to high z and possess scaling relations with as small 
as possible intrinsic scatter, but its physics should be well 
understood from a theoretical point of view. On the con- 
trary, there is up to now not any definitive understanding 
of the GRBs scaling relations so that, as a dangerous draw- 
back, one can not anticipate whether the calibration param- 
eters are redshift dependent or not. Lacking any theoretical 
model, we will therefore address here this problem in a phe- 
nomenological way adopting two different parameterizations 
for their evolution with z allowing for a change in the zero- 
point only or in both zeropoint and slope. It is worth noting 
that such an analysis has been severely hampered before by 
the low statistics in the GRBs sample, a problem that we 
partially overcome here thanks to the use of the 115 GRBs 
catalog recently assembled by Xiao & Schaefer (2010). 

The plan of the paper is as follows. In Sect. 2, we briefiy 
review the different GRBs 2D correlations known insofar and 
present the Bayesian motivated method we will use in the 



following to calibrate them using three different approaches 
to estimate the GRBs luminosity distances. The constraints 
thus obtained and their dependence on the adopted luminos- 
ity distance determination are discussed in Sect. 3, while, in 
Sect. 4, we address the problem of their dependence on the 
redshift. The issues related to the HD derivation are then 
discussed in Sect. 5 where we consider the impact on the 
HD of the luminosity distance estimate method, of the the 
averaging procedure and the homogeneity of the sample. A 
summary of the main results and a discussion of them is 
finally given in Sect. 6, while some complementary material 
is presented in Appendix. 



2 GRBS SCALING RELATIONS 

It is instructive to start this analysis considering the general 
case of two observable quantities {x, y) related by a power - 
law relation which, in a log -log plane, reads 



log y = a log X + b 



(1) 



Calibrating such a relation means determining the slope a, 
the zeropoint b and the scatter ai„t of the points around 
the best fit relation. In order to be useful for distance de- 
termination, one of the two quantities, say x, should refer 
to a directly observed quantity, while the other one must 
depend on the redshift z. Setting y = Kdj^{z) with k a di- 
rectly measurable redshift independent quantity and (Il^z) 
the luminosity distance, we get 

log y = log K -f 2 log cLl {z) = a log a; -I- 6 

so that one can then estimate the distance modulus as : 

= 25 + ^\ogdL{z) = 25 + (5/2) (aloga; + & - logs) . 

In order to perform such an estimate, a two steps procedure 
has to be implemented. First, one has to select a sample 
of low redshift (z ^ 0.01) objects with known distance and 
fit the scaling relation to the [x, y) data thus estimated to 
infer the calibration parameters {a,b,aint)- Second, one has 
to assume that such calibration parameters do not change 
with the redshift so that a measurement of (x, k, z) and the 
use of the above scaling relation, with {a,b,aint) the same 
at all z, are sufficient to infer the distance modulus. 

Both these steps are daunting tasks for GRBs because 
of the lack of low redshift objects to be used in the calibra- 
tion procedure. In order to overcome this problem, different 
strategies can be implemented, but their impact on the fi- 
nal derivation of the GRBs Hubble diagram have not been 
investigated in detail which is our aim here. 



2.1 2D empirical correlations 

Their high luminosity offering the possibility to be detected 
at very large z makes GRBs promising candidates to trace 
the Hubble diagram in the matter dominated era so that a 
great effort has been devoted during the late years to look 
for reliable and narrow empirical correlations. We limit here 
our attention only to two dimensional (hereafter, 2D) cor- 
relations since they can be investigated relying on a larger 
number of GRBs. These involve a wide range of GRBs prop- 
erties related to both the energy spectrum and the light 
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curve which are then correlated with the isotropic luminos- 
ity L or the emitted coUimation corrected energy E.^. Nei- 
ther L nor are directly measurable quantities since they 
depend on the luminosity distance dL^z). Indeed, it is: 



E^f = ^■Kd'i(z) SboloFbeam.il + z) 



(2) 
(3) 



where Pboio and Sboio are the bolometric peak flux and flu- 
ence, respectively, while Fbeam = 1— cos (Ojet) is the beaming 
factor with 6 jet the rest frame time of achromatic break in 
the afterglow light curve. Note that the bolometric qua nti- 
ties are related to the observed ones as (|Schaefer II2OO7I ) : 



P X 



S X 



W/Ci+z) 

/(1+^) 



E${E)dE 



J^""- E${E)dE 



(4) 



(5) 



with P and S the observed peak energy and fluence and 



E max 



the detection thresholds of the observing 



instrument. We model the ener gy spectrum ^(E ) as a 
smoothly broken power - law, i.e. (|Band et al. II 1993 1 : 



$(£;) 



AE° 



BE'' 



exp 



{2 + a)E 



E, 



peak 



E 



E; 



peak 



2 + a 



(6) 



otherwise 



The GRBs 2D correlations are obtained setting y — L or 
y = E.y in Eq.([l]), while different choices are available for x. 
In particular, we will consider the following possibilities : (i.) 
X = Epeak{l + z)/300 wlth Epeak the peak energy (in keV) 
of the fluence spectrum; (ii.) x = riagil + 2)^^/0.1 with riag 
(in s) the time offset between the arrival of the low and high 
energy photons; (iii.) x = r_RT(l + z)~^ /O.l with trt (in 
s) the shortest time over which the light curve rises by half 
the peak flux of the burst; (iv.) x = V(l + z)/0.02 with V 
the variability which quantifles the smoothness of the light 
curve itself. Note that all these quantities are expressed in 
the GRB rest frame which motivates the redshift term, while 
the further scaling constant is introduced to minimize the 
correlation among errors. 

The combination of x and y gives rise to 
different co rrelation we w il l con s idefl namely 
E.y - E„e nk (ichirlanda et alTI 12004 I2OO6I fh . the 



Epenh (|Schaefer]~ 2OO3I : Yonekotu et al 



the 
the 
L- 



200d) 



^, 200i l. L-Tlag 

L-TRT JSchaefer I |2007|) and L-V 



Norris et al. . _ , .. . . 

I Fenimore fc Ramirez - Ruiz I I2OO0I : iRiechart et al.1 I2OOII : 



ISchaeferTl2007l '). These flve correlations has been used in 
Schaefer (2007, hereafter S07) to derive a combined Hub- 
ble diagram for 69 GRBs, later updated by Cardone et al. 
(2009, hereafter CCD09) where also a sixth correlation be- 
tween the break time Ta and t he luminosity at break time 
Lx(T a) of the X-ray afterglow (jPainotti et al. II2OO8I . [ioiol . 
I2OIII ) has been considered. 

Recently, a larger GRBs catalog has been compiled by 
Xiao & Schaefer (2009, hereafter XS09) giving, for each 



^ Hereafter, we will refer to the correlation obtained setting x = 
X and y = Y in Eq.l[T]| as the Y - X correlation or scaling law. 



object, the values of the different quantities (if available) 
needed to check the flve correlations above. We use this 
sample here as input for our analysis referring the inter- 
ested reader to XSIO for details on the catalog construction 
and observable quantities determination. 

2.2 Bayesian fitting procedure 

Eq.Q is a linear relation which can be fltted to a given 
dataset (xi,yi) in order to determine the two calibration 
parameters {a,b). Moreover, although there is still not a 
theoretical model explaining any of the empirical 2D cor- 
relations in terms of GRB physics, we nevertheless expect 
that the wide range of GRB properties makes the objects 
scatter around this (unknown) idealized model. As a con- 
sequence, the above linear relations will be affected by an 
intrinsic scatter (Jint which has to be determined together 
with the calibration coefficients {a,b). To this aim, in the 
following we will r esort to a Bayesian motivated technique 
ID' Agostinil (|2005l ) thus maximizing the likelihood function 
£(a, b, Oint) — exp [—L{a, b, aint)] with : 



L(a, b, aint) 



1 (Yi- aXi - bf 



j2 

ZTlt 



4, 



(7) 



with {Xi,Yi) = {\ogXi,\ogyi) and the sum is over the A/" 
objects in the sample. Note that, actually, this maximization 
is performed in the two parameter space (a, aint) since b may 
be estimated analytically as : 



+ ai., + a^ai 



+ ai-, + a^ai 



(8) 



so that we will not consider it anymore as a fit parameter. 

It is worth noting that the Bayesian approach allows to 
flnd out what is the most likely set of parameters within a 
given theory, but does not tell us whether this model flts 
well or not the data. An easy way to quantitatively estimate 
the goodness of the flt is obtained considering the median 
and root mean square of the best fit residuals, defined as 
5 = Yobs — Yfit which we will also compute for the different 
2D correlations we will consider. 

The Bayesian approach used here also allows us to quan- 
tify the uncertainties on the fit parameters. To this aim, for 
a given parameter pi , we first compute the marginalized like- 
lihood Ci{pi) by integrating over the other parameter. The 
median value for the parameter pi is then found by solving : 



Ci{pi)dpi 



Ci{pi)dpi 



(9) 



The 68% (95%) confidence range {pi^i,Pi^h) are then found 
by solving : 



Ci{pi)dpi 



Ci{pi)dpi 



Ci{pi)dpi 



Ci{pi)dpi 



(10) 



(11) 



with e = 0.68 (0.95) for the 68% (95%) range respectively. 
Actually, in order to sample the parameter space, we use a 
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Markov Chain Monte Carlo (MCMC) method running two 
paraUel chains and using the Gelman - Rubin (1992) test 
to check convergence. The confidence ranges are then ob- 
tained considering the histograms of the parameters from 
the merged chain after burn in cut and thinning. 



2.3 GRBs luminosity distances 

A preliminary step in the analysis of all the 2D correlations 
hinted at above is the determination of the luminosity L 
or the coUimated energy E-y entering as Y variable in the 
Y-X scaling laws. As shown by Eqs.Q - Q, one has to 
first determine the GRBs luminosity distance over a redshift 
range where the linear Hubble law does not hold anymore. 
As such, one should estimate the luminosity distance as : 



dL{z) 



C{1 



Ho 



dz' 



(12) 



with Ho = lOOh km/s/Mpc the present day Hubble constant 
and E{z) = H{z)/Ho the dimensionless Hubble parameter 
which depends on the adopted cosmological model thus lead- 
ing to the well known circularity problem (i.e., one would like 
to use GRBs to probe the cosmological model, but needs the 
cosmological model itself to get the GRBs Hubble diagram). 

Different strategies have been developed to tackle this 
problem. The simplest one is to assume a fiducial cosmolog- 
ical model and determine its parameters by fitting, e.g., the 
SNela Hubble diagram. Motivated by its agreement with a 
wide set of data, one usually adopt the ACDM model as 
fiducial one thus setting : 



E\z) = nM{i + zf 



(13) 



with SIa = I — Q,M because of the spatial fiatness assump- 
tion. In order to determine the parameters {Qm, h), we max- 
imize the likelihood function £,{Qm, h) oc exp (— x'^/2) with : 



E 



+ 



h°^' - h 
CTh 



(14) 



As input data, we use the Union2 SNela sample 
lAmanuUah et aL I [ioiol ') to get (^f^cr^,) for NsNeia = 
557 objects over the redshift range (0.015, 1.4) and set 
(a;M^o■^A/) = (0.1356,0.0034) for the matter physical den- 
sity UJM = ^Mh^ and {h,(Jh) = (0.742,0.036) for the 
Hubble c onstant in agreement with the re sults of the 
WMA P7 jKomatsu et al. II2OIII ) and SHOES (|Riess et al. I 
I2OO9I ) teams, respectively. The best fit values turn out to be 
{Qm, h) = (0.261, 0.722), while median value and confidence 
ranges read : 

O — n OKQ+0 019 +0-032 , 
\Lm — U./t)y_o.oi6 -0.028 1 " 



f^ 700+0.025 +0.046 
U. ( ^O_0.025 -0.045 



Although the ACDM model fits remarkably well the data, 
it is nevertheless worth stressing that a different cosmolog- 
ical model would give different values for dL^z) and hence 
different values for {L,Ej) thus impacting the estimate of 
the calibration parameters {a,crint). Although CCD09 have 



shown that this effect is quite small for a large range of phe- 
nomenological dark energy model^, it is nevertheless worth 
be conservative and look for model independent approaches. 
A first step towards this aim is to resort to cosmography 
j Weinberg II 1972l :lv isser |[20o3 ) . i.e., to expand the scale fac- 
tor a{t) to the fifth order and then consider the luminosity 
distance as a function of the cosmographic parameters. In- 
deed, such a kinematic approach only relies on the validity 
of assumption of the Robertson - Walker metric, while no as- 
sumption on either the cosmological model or the theory of 
gravity is needed since the Friedmann equations are never 
used. We use the expansion of dL{z) to the fifth order in z 
found in Capozziello et al. (2008) and determine the Hubble 
constant h, the deceleration qo, jerk jo, snap so and lerk lo 
parameters by fitting the Union2 SNela data with the prior 
on h coming from the SHOES team quoted above. Note that, 
in this case, the likelihood function is defined as before, but 
the pseudo-x^ is now given by Ea. p4|) without the term 
depending on ujj\i since now is no more a parameter. 
Using a similar MCMC algorithm (but now running three 
parallel chains) , we find as best fit values : 

{h, go, jo, so, lo) = (0.741, -0.56, 0.66, -0.41, 3.59) , 



while median values and confidence ranges read : 

10 



h — n 7/11+0.036 +0.071 
n — U.(y:i_o.o35 -0.072 



^(-+0.05 +0.10 
'U.*0_o.05 -0.14 



n nn+013 

JO = O.OOI0.12 



-0.79 
-0.35 



a — n OQ+0-74 +1.57 
So — U.Zy_o.20 -1.08 



-0.07 



+4.77 +12.27 



in good agreement wit h prev io us results i n liter- 
ature (IVitagliano etlTI I2OI0I : IXu fc Wangi I2OI0I : 
iBouhmadi - Lopez et al. II2OIOI ) using different datasets. 

As a further step towards a fully model independent es- 
timate of the GRBs luminosity distances, one can use SNela 
as distance indicator based on the naive observations that a 
GRBs at redshift z must have the same distance modulus of 
SNela having the same redshift. Interpolating therefore the 
SNela Hubble diagram gives therefore the value of jj,{z) for a 
subset of the GRBs sample with z ^ 1.4 which can then be 
used to calibrate the 2D correlations ( Kodama et al. I I2OO8I : 
iLiang et al."1l2008l : I Wei fc Zhang 1120091 ). Assuming that this 
calibration is redshift independent, one can then build up 
the Hubble diagram at higher redshifts using the calibrated 
correlations for the remaining GRBs in the sample. CCD09 
have used a si milar approach based on the local regres- 
sion techniqu e (|Cleveland II 19791 : ICleveland fc Devlin 1 ll988l : 
iLoader ll 19991 ) which combines much of the simplicity of lin- 
ear least squares regression with the flexibility of nonlinear 
regression. The basic idea relies on fitting simple models to 
localized subsets of the data to build up a function that de- 
scribes the deterministic part of the variation in the data, 
point by point. Actually, one is not required to specify a 



2_CCD09 adopted the CPL i jChevallier fc Polarski ll200ll:lLinder1 
l2003t) parametrization of the dark energy equation of state, w = 
wo + Waz/{1 + 2), and explored the dependence of {a,crint) on 
{wq, Wa) for different values of Qm- While clear trends of (a, aint) 
with {Qm,wo,Wa) were obtained, the relative change is negligi- 
bly small. It is worth stressing, however, that the calibration pa- 
rameters may change more if the adopted cosmological model is 
radically different from the ACDM one (sec, eg, Diafcrio et al. 
2011 for an illuminating example). 
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global function of any form to fit a model to tfie data so 
that there is no ambiguity in the choice of the interpolating 
function. Indeed, at each point, a low degree polynomial is 
fit to a subset of the data containing only those points which 
are nearest to the point whose response is being estimated. 
The polynomial is fit using weighted least squares with a 
weight function which quickly decreases with the distance 
from the point where the model has to be recovered. We ap- 
ply local regression to estimate the distance modulus 
from the Union2 SNela sample following the steps schemat- 
ically sketched in CCD09 which we refer the reader to for 
further details and the demonstration of the reliability of 
the inferred luminosity distances. 



3 CALIBRATION PARAMETERS 

While the X quantities are directly observed for each GRBs, 
the determination of Y (either the luminosity L or the col- 
limated energy E.y) needs for the object's luminosity dis- 
tance. The three methods described above allows us to get 
three different values for Y so that it is worth investigating 
whether this has any significant impact on the calibration 
parameters {a,b,ai„t) for the correlations of interest. We 
will refer hereafter to the three samples with the Y quanti- 
ties estimated using the luminosity distance from the fidu- 
cial ACDM cosmological model, the cosmographic parame- 
ters and the local regression method as the F, C and LR 
samples, respectively. 

As a preliminary caveat, it is worth stressing that the 
error on Y comes from two different contributions. First, 
there is the uncertainty obtained by propagating those on 
the observed quantities (i.e., the fiux P and the Band pa- 
rameters for L and the fluence S, the jet angle 9 jet and again 
the Band parameters for E^,) which we assume to be uncor- 
related. Second, there is the uncertainty on dL{z) which we 
estimate from the MCMC chain for the F and C cases, while 
it is output from the method for the LR case. It is worth 
noting that, given the large SNela sample used in the fit of 
ACDM and cosmographic parameters and in the local re- 
gression technique, the error on the luminosity distance is 
actually quite small and always smaller than the one from 
the measurement uncertainties. It is this latter term that 
one should minimize in order to get better determined L 
and/or E.y values and hence put stronger constraints on the 
slope, zeropoint and intrinsic scatter. 

In order to constrain the calibration parameters 
{a, h, (Tint), we can use the Bayesian procedure described 
above with the three GRBs samples as input to the like- 
lihood analysis. However, the F, C and LR samples can 
not contain the same number of objects since the GRBs 
luminosity distance (and hence the Y quantities) can be es- 
timated only for Zmin ^ 2 ^ Zmax with {zmimZmax) de- 
pending on the adopted method. If we rely on the fidu- 
cial ACDM model, we can predict dL{z) at every z so 
that {zmin, Zmax) = (0, oo) with the only caveat that one 
is implicitly assuming that a model fitted over the range 
(0.015, 1.4) probed by the Union2 SNela sample can be ex- 
trapolated to the full evolutionary history of the universe. 
The C sample is based on the use of the cosmographic pa- 
rameters fitted to the data having checked that the fifth 
order expansion of dL{z) is reliable over the Union2 SNela 



redshift range. However, the larger is z, the higher is the 
order one has to include in the dL{z) Taylor series to get 
accurate results so that one must set {zmin, Zmax) ~ (0, 1.4) 
in order to not bias the determination of the Y quanti- 
ties because of an inaccurate distance approximation. Fi- 
nally, the local regression method can only be applied to 
estimate dL{z) over the redshift range probed by the data 
used for the interpolation so that one get the constraint 
{zmin, Zmax) ~ (0.015, 1.4). lu Order to homogenize the three 
samples, we therefore fit the 2D correlations using only F, 
C and LR GRBs with 2 1.4 so that we can meaningfully 
compare the results for the three distance estimate methods. 

Table[T] summarizes the constraints on the calibration 
parameters for the three different samples. Note that we 
have also included the Eiso - Ep^ak correlation (|Amati et al. I 
l2009t ) with E,so = E^/Fieam since, although refers to the 
same quantities used in the E-^ - Epeak correlation, it is based 
on a larger number of GRBs (given that one does not need 
a measurement of the jet opening angle 6 jet)- As a general 
result, we find that the fit is always quite good with reduced 
values close to 1 in all the cases independently on the 
2D correlation considered and the distance estimate method 
adopted. As such, the results from the different sample are 
statistically equivalent and can be safely compared. 

Although the 68 and 95% confidence ranges for the cal- 
ibration parameters {a,b,aint) for a given correlation over- 
lap quite well for the F, C and LR samples, the best fit 
coefficients and the median values clearly show that the cal- 
ibration based on the fiducial ACDM model leads to steeper 
scaling laws for most of the cases. On the contrary, shal- 
lower slopes are obtained using the C or LR samples with 
the L-V relation as unique exceptioijf]. Although the differ- 
ences in the slopes are not statistically meaningful because 
of the large uncertainties, it is nevertheless worth investigat- 
ing whether such an effect can be ascribed to the distance 
determination. To this end, for a given GRB and Y -X cor- 
relation, let us define the following quantities : 

AxY = (aiX + h) - {a^X + fea) = (ai - a2)X + (bi - fea) , 

Aa^Y ^2log[dLl{z)/dL2{z)] . 

While AxY gives the difference in the Y values predicted 
using the correlations with coefficients (ai,6i) and (a2,&2), 
Adj^Y quantifies the effect of estimating the model depen- 
dent quantity Y using two different luminosity distances. 
Should the change in the slope be the outcome of compen- 
sating the offset due to different luminosity distances as- 
sumptions, one should get AxY ~ Adj^Y. Actually, this is 
not the case. Considering, e.g, the L-Ep^ak correlation, a 
weighted average gives : 

(AxF) =0.19/ (Ad^y) = 0.09 

using the F and C samples and 

(AxY) = -0.25 / (Ad^Y) = -0.02 

for the F vs LR samples. We therefore find that the change 

^ Actually, the L-V relation is the shallowest one and has the 
largest intrinsic scatter among the six 2D correlations investigated 
in this paper so that a different trend of the slope with the sample 
could also be a statistical artifact. 
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Table 1. Constraints on the calibration parameters (a,b,aint) for the 2D correlations considered in the text. Columns arc as follows: 
1. correlation id, 2. number of GRBs used, 3. reduced x^i 4. best fit parameters, 5., 6., 7. median values and 68% confidence ranges for 
(a, 6, (Tint), respectively. For each correlation, there are three rows referring to the results obtained using the F, C and LR distances. 



Id 



M xVd.o.f. 



(a, b, (7int)bf 



L-E, 



peak 



L-E, 



peak 



L-E, 



'peak 



39 
39 
39 



1.05 
1.04 
1.08 



1.05, 49.75, 0.46) 

1.06, 49.66, 0.47) 
0.95, 50.00, 0.43) 



0.99 
1.02 
0.84 



-0.21 +0.41 
-0.22 -0.43 

+0.21 +0.44 
-0.20 -0.43 

+0.23 +0.47 
-0.23 -0.44 



49.61 



hO.64 +1.14 
-0.21 -0.97 



0.46"* 



hO.16 
-0.11 



49.53 



+0.57 +1.15 



0.47 



+0.07 +0.15 



50.02^ 



1-1.24 
-0.72 



0.45 



+0.05 +0.17 



L-Tiag 27 1.07 

L-r,,g 27 1.09 
L-Ti,„ 27 1.09 



0.70, 51.60, 0.47) 
0.68, 51.57, 0.47) 
0.67, 51.67, 0.43) 



-0.65 
-0.64 
-0.62 



+0.17 +0.36 
0.16 -0.31 

+0.14 +0.28 
-0.15 -0.33 

+0.15 +0.27 
0.14 -0.33 



51.60 
51.60 
51.70 



+0.12 +0.25 
0.05 -0.17 

+0.07 +0.17 
-0.08 -0.21 

+0.10 +0.18 
0.07 -0.18 



0.48 

0.49 
0.44 



+0.06 +0.19 
0.06 -0.11 

+0.08 +0.20 
-0.07 -0.13 

+0.07 +0.16 
0.06 -0.12 



L-TRT 31 1.16 

L-TRT 31 1.20 
L-tut 31 1.10 



1.09, 51.81, 0.45) 
1.03, 51.78, 0.46) 
1.07, 51.88, 0.40) 



-1.02 



+0.24 +0.46 
-0.19 -0.47 



-0.92 



+0.25 +0.47 
-0.24 -0.49 



-0.97 



+0.22 +0.45 
-0.22 -0.47 



51.82 
51.83 
51 



+0.05 +0.12 
0.04 -0.10 

+0.04 +0.10 
-0.07 -0.15 

+0.07 +0.14 
-0.02 -0.10 



0.46 
0.48 
0.41 



+0.08 +0.20 
0.06 -0.11 

+0.08 +0.18 
-0.07 -0.13 

+0.07 +0.15 



L-V 
L-V 
L-V 



37 
37 
37 



1.10 
1.09 
1.10 



0.34, 52.60, 0.61) 
0.40, 52.65, 0.61) 
0.37, 52.73, 0.56) 



0.33 



+0.18 +0.43 
-0.20 -0.52 



er, co + 0.25 +0.56 r, Rr, + 0.10 +0.20 

0^-00_0.16 -0.47 ^■°^-0.07 -0.14 



0.37 



0.21 +0.49 
-0.25 -0.50 



52.66 



+0.16 +0.54 
0.30 -0.69 



0.63 



+0.09 +0.23 
0.08 -0.14 



0.33 



-0.19 +0.41 
-0.22 -0.52 



7+0.09 +0.13 



Ej-Ep^ak 14 1.17 



Ej - Ep^ak 



14 



E~f - Ep^ak 14 



1.21 
1.16 



1.77, 46.71, 0.17) 

1.71, 46.80, 0.18) 

1.72, 46.81, 0.17) 



1.71 



+0.22 +0.44 
-0.24 -0.65 



,+0.23 +1.01 



-,+0.07 +0.19 



1.62 



+0.24 +0.51 
-0.29 -0.64 



47.28 



+0.27 +1.10 
-0.71 -1.35 



0.21 



+0.07 +0.16 



1.66" 



■0.28 -0.57 



q+0.3S +0.95 



-1+0.07 +0.18 



Eiso-Ep^ak 40 0.79 
Eiao - Ep^ak 40 0.81 

Eiso-Ep^ak 40 0.79 



1.58, 49.16, 0.52) 
1.24, 49.57, 0.52) 
1.22, 50.03, 0.51) 



1.54 



+0.27 +0.55 
-0.26 -0.50 



49.26 



+0.45 +1.03 
-0.45 -1.04 



0.51 



+0.09 +0.21 



1.00 



+0.37 +0.75 
-0.33 -0.62 



50.16 



+0.99 +1.69 
-0.30 -1.27 



0.52 



+0.11 +0.24 



1.04 



+0.35 +0.73 
-0.38 -0.70 



3+0.97 +1.77 



,+0.12 +0.27 



in the slope is not induced by the different luminosity dis- 
tances adopted. Actually, the use of different samples has 
also an impact on the intrinsic scatter determination with 
the LR sample leading to smaller best fit and median aint 
values for all the six correlations considered. Since a and 
aint correlate, the change in the slope is not only due to the 
change in the luminosity distances, but also of the intrinsic 
scatter. Discriminating among these two effects and quanti- 
fying their respective contribution is not possible so that we 



can not draw any definitive conclusion on the impact of the 
calibration method on the slope of the 2D scaling relations. 



4 EVOLUTION WITH REDSHIFT 

As more and more data add to the available GRBs dataset, 
the observational evidences for the six 2D correlations we are 
considering become more and more reliable. On the contrary, 
a clear theoretical motivation is still lacking in many cases. 
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Table 2. Constraints on the calibration parameters (a,b,a,cri„t) for the 2D correlations considered in the text and fitted to the F 
sample. Columns arc as follows: 1. correlation id, 2. reduced x'^i 3. best fit parameters, 4., 5., 6., 7. median values and 68% confidence 
ranges for the fit parameters. The number of GRBs used in the same as in Table[T] 



Id x^/d-o-f- {a,b,a,aint)bf aba ai„t 





1.08 


(0.89, 50.77, 0.67, 0.41) 


0.82±»;?° 


50, 


,75; 


f 0,59 
-0.65 


p, fio+0,18 




L-Tlag 


1.10 


(-0.58, 50.38, 0.88, 0.39) 


-0.57t«;- 


50, 


,39; 


f0,66 
-0.58 


n. 07+0,30 


^^•^^-0,06 


L-TUT 


1.16 


(-0.79, 51.22, 0.65, 0.44) 




50, 


,97; 


f0,70 
-0.75 


'"-0.33 


47+0.07 

-0,06 


L-V 


1.23 


(-0.10, 50.27, 1.07, 0.54) 


-0.05±«;f, 


50, 


,29; 


f0,91 
-0.96 


1 04+'''^2 


'-'•'"-0,08 


- Ep^a}^ 


1.37 


(1.57, 50.64, 0.00, 0.14) 


1 48+°'^^ 


50, 


42+0,51 


0.131^2^ 


"•^'J-0,06 


^iso ~ Epeak 


0.81 


(1.39, 51.58, 0.58, 0.48) 


1 oq+O^S 
-^-•^"-0.28 


51, 


,57; 


f0,84 
-0.81 


0.581^4^ 


'J-O'J-0,08 



Table 3. Constraints on the calibration parameters (ag, log ai , boi log ''1 1 f in* ) for the 2D correlations considered in the text and fitted 
to the F sample. Columns arc as follows: 1. correlation id, 2. reduced x^i 3. best fit parameters, 4., 5., 6., 7., 8. median values and 68% 
confidence ranges for the fit parameters. The number of GRBs used in the same as in Tablc[l] 



Id x^/d.o.f. (ao, log ai, 60, log 61, o-i„t)6 J ao logai 60 logfei (Tint 



^ Epeak 


1.12 


(0.89, -2.69, 51.49, -0.21, 0.41) 


0-89to,25 


'^■^''-3.63 


51 R7+0-^S 
Cil-«'_0,46 




44+0-03 
"•^^-0.05 


L-Tiag 


1.13 


(-0.60, -4.43, 51.31, -0.07, 0.40) 


-0 64+°-l'^ 


70 + 1,61 
-^■'»_2,85 


51 fil+0-49 
ai.Oi_g 45 


-0.28l«j^ 


43+0-08 
^•^■'-0.06 


L-TRT 


1.17 


(-1.56, -0.08, 52.38, -0.70, 0.44) 


-0 95+°-^3 

'^■^''-0,46 


1 0-1+1,59 
-1-»1_3.90 




1 SI +1-60 


47+0-07 

'^•^'-0.06 


L-V 


1.36 


(-0.15, -2.16, 51.40, 0.03, 0.52) 


-0 10+°-*^ 


-l-79ti-i? 




_o 77+0,76 
'^•' '-2,86 


n co+o.io 

0-&»_o.08 


E^ - Ep^Q^f^ 


1.70 


(0.79, 0.04, 50.37, -0.56, 0.10) 


-'-■^'-'-0,47 


-1.70l^;f7 


50.6ltr^ 


-2.86j;3;g4 


Q 1 7+0.08 
'^•l'-0.06 


^iso ~ ^peak 


0.79 


(0.94, -0.23, 52.00, -0.15, 0.49) 


-I 97+0,33 
l-^'-0,50 


-1.78l-«l 


52.56l«-« 


-I 70+1,42 
l-'-'-2.41 


^•"^-0.08 



As a consequence, it is also not clear whether the calibration 
parameters {a,b,aint) evolve with the redshift or not. To 
investigate this issue, we consider two different possibilities 
for the evolution with z. First, we consider the possibility 
that the slope is constant, but the zeropoint is evolving. In 
particular, we assume : 

y = B{1 + z)"x'^ — > Y = alog{l + z)+aX + b (15) 

with (X,Y) = (log a;, logy) and {a,b) = (A,log&). Tableg] 
summarizes the constraints on the {a,b, a, ai„t) parameters 
obtained fitting Eq. (|15|) to the 2; ^ 1.4 GRBs with distances 
estimated from the fiducial ACDM model. The constraints 
from the fit to the C and LR samples are consistent within 
the 68% confidence ranges so that they are reported in Ap- 
pendix for completeness, but not discussed anymore here. 

Comparing to the constraints in Table[l]highlights some 
interesting lessons. First, we note that both the best fit and 
median values of the slope parameter a are significantly shal- 
lower than in the no evolution case. However, the 68% con- 



fidence ranges typically overlap quite well so that, from a 
statistical point of view, such a result should not be over- 
rated. Actually, adding one more parameter introduces a 
degeneracy between a and a so that one can make the best 
fit K-X relation shallower compensating the difference in 
the term aX with the contribute from a{l + z). Indeed, a is 
typically quite large the only exception being the E-y - Epeak 
correlation which is also the only one with the same best fit 
slope for the fit with and without the evolution term. Quite 
surprisingly, the no evolution result (i.e., a = 0) is consis- 
tent with the 68% confidence range only for the E-y - Epeak 
correlation thus arguing in favor of an evolution of the cali- 
bration parameters with 2. On the other hand, the reduced 
values are never smaller than those obtained for the fits 
with no evolution. Moreover, the intrinsic scatter is almost 
the same for both cases so that allowing for a redshift evolu- 
tion does not lead to statistically preferred results. As such, 
we consider a most conservative option to assume that the 
GRBs scaling relations explored here do not evolve with z. 
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Actually, such a conclusion is model dependent. As an 
alternative parametrization, we therefore allow for an evo- 
lution of the slope and not only the zeropoint of the 2D 
correlations. We fit the data using : 



y = (ao + aiz)X + {bo + hiz) 



(16) 



i.e., we are Taylor expanding to first order the unknown 
dependence of the slope and zeropoint on the redshift. Note 
that, actually, we should limit the redshift range to very 
low z to have a meaningful expansion, but we extrapolate 
this linear relation to any z as a first order guess avoiding 
to add further unknown fitting parameters. Morevoer, since 
we expect (ai,foi) are quite small, we skip to logarithmic 
units for these quantities in order to be more sensitive to 
the tiny values. Table[3] summarizes the results for the fit to 
the F sample, while a qualitatively similar conclusions can 
be drawn from the those to the C and LR samples as can 
be inferred from Table[2]in Appendix A. 

As a general result, we find that the best fit parame- 
ters and the median values of the evolutionary coefficients 
(log oi, log 61) are typically quite small indicating that the 
dependence of both the slope and the zeropoint on the red- 
shift is quite weak, if present at alQ . On the other hand, 
the slope and zeropoint at 2; = are consistent within the 
68% confidence range with the corresponding quantities in 
the no evolution case. Moreover, the reduced values are 
typically larger than in the no evolution case so that this 
latter assumption is statistically preferred. 



5 GRBS HUBBLE DIAGRAM 

Once the calibration parameters for a given Y -X correla- 
tion have been obtained, it is then possible to estimate the 
distance modulus of a given GRB from the measured value 
of X. Indeed, for a given Y , the luminosity distance is : 



dl{z) 



with K = A-rrPboio, « = 4:TvSboioFbeam/{'i- + z) and k = 
47rSi,oio/(l + z) for Y = L, Y = Ej and Eiso, respectively. 
Using the definition of distance modulus and estimating Y 
from X through the Y -X correlation, we then get : 



fi{z) = 25 + 51ogdi(2) 

= 25 + (5/2) log (y^) 

= 25 + (5/2)(aA: + &-logK) 



(17) 



where {a,b) are the best fit coefficients for the given Y-X 
correlation. Note that we will refer here only to the fits with 
no redshift evolution of the calibration parameters since the 
analysis in the previous section has not shown any clear ev- 
idence for a redshift dependence of {a,b). Although such a 
result is model dependent (since we have considered only 



* Note that, since we are using logarithmic units, the no evo- 
lution case can not be exactly achieved corresponding to log ai 
and logfei going to —00. However, needless to say, obtaining, e.g., 
logai ~ —2 is the same as stating that there is no evolution at 
all. Note also that we have implicitly assumed that both ai and 
61 are positive. We have, however, checked that the qualitative 
conclusions are not affected by this assumption. 



two possible evolution parameterization) and actually lim- 
ited only to the redshift range 0.015 ^ 2 ^ 1.4, we are nev- 
ertheless confident that the possible bias induced by the no 
evolution assumption is smaller than the present day mea- 
surement uncertainties. However, such an issue worths to be 
reconsidered when a larger and more precise sample will be 
available. 

Ea. (|17|l allows to compute the central value of the dis- 
tance modulus relying on the measured values of the GRBs 
observables, i.e., the ones entering the quantity k, and the 
best fit coefficients (a, b) of the used correlation. However, 
both K and (a, b) are known within their own uncertain- 
ties which have to be propagated to get the error on fJ.(z). 
Moreover, each correlation is affected by an intrinsic scat- 
ter which has to be taken into account into the total error 
budget. To this end, we adopt the procedure schematically 
sketched below. 

(i) Fix (a, (Tint) to the i-th value of the Markov Chain 
and estimate b from Eq.(|S}. 

(ii) Set the quantities needed to compute X and hence k 
randomly sampling from Gaussian distributions centred on 
the measured observed values and with variance equal to the 
measurement error. 

(iii) Compute fi from Ea. (fT7)) for aU the 2000 X values 
generated above and estimate the mean of the distribution. 

(iv) Repeat steps (ii) and {iii) for all the points along the 
chain and estimate the statistical error on fi by symmetrizing 
the 68% confidence range of the distribution thus obtained. 

(v) Add in quadrature the statistical error and the in- 
trinsic scatter ai„t of the correlation to finally get the total 
uncertainty on the distance modulus fj,. 

Eq. (|17|l and the above procedure allows us to build up the 
Hubble diagram for all the GRBs with measured values of 
the quantities entering a given Y-X correlation. It is then 
possible to both reduce the uncertainties and (partially) 
wash out the hidden systematic errors by averaging over the 
different correlations available for a given GRB. Following 
Schaefer (2007), we finally estimate the distance modulus 
for the i - th GRB in the sample at redshift Zi as : 



E 



l/c 



with the uncertainty given by : 



(18) 



(19) 



where the sum runs over the 2D empirical laws which can 
be used for the GRB considered. 

As summarized in Table [1] and discussed in Sect. 3, the 
best fit calibration parameters depend on the method used 
to set the GRBs luminosity distances. Although such a de- 
pendence is not statistically meaningful being the confidence 
ranges for the different cases well overlapped, it is neverthe- 
less worth investigating whether they have an impact on 
the derived Hubble diagram. Moreover, averaging over more 
than one correlation implicitly assumes that all of them are 
physically motivated (i.e., they are the outcome of an un- 
known underlying theoretical mechanism) and not the arti- 
fact of some not well understood selection effect. It is there- 
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Figure 1. GRBs Hubble diagrams (HDs) averaging over the six 2D correlations and exeluding outliers (see text for the definition 
of outliers). Three panels refer to the HDs derived using the calibration based on the fiducial ACDM model (left), the cosmographic 
parameters (centre), and the local regression (right). 



fore also important to check what should be the effect of 
a possible misleading assumption. Both these issues will be 
discussed in the following. 

5.1 Impact of the calibration method 

Fig.[T] shows the GRBs Hubble diagrams (hereafter, HDs) 
obtained averaging over the six 2D correlations and using 
the three different calibration methods. As a guidance for 
the eye, the red solid line is the expected /^(z) curve for the 
fiducial ACDM model. Note that we have excluded GRBs 
which deviate from this line more than la a with a a the rms 
of fifid — fJ.{z) for the full sample. Such a criterium is actually 
quite loose excluding only two objects out of 114 so that we 
are confident that no bias is induced by this cut. 

As a general remark, we find that, notwithstanding the 
calibration method adopted, the GRBs HDs reasonably fol- 
low the ACDM curve although with a non negligible scatter. 
Quite surprisingly, the scatter is significantly larger in the 
range 0.4 ^5 z ^ 1.4 because of a set of GRBs with fi{z) ly- 
ing systematically above the ACDM prediction. One should 
argue for a failure of the theoretical model, but there are 
actually a set of points which are hard to reconcile with any 
reasonable dark energy model. While this could be dropped 
out from the sample by an iterative selection procedure, we 
prefer to err on the conservative side in order to introduce 
any bias induced by an a priori choice of a dark energy 
model. It is nevertheless likely that these GRBs are outliers 
of one or more of the 2D scaling laws used to infer their 
distance modulus. Looking in detail to their ij,{z) estimates, 
we find that the estimates obtained from the different corre- 
lations entering the averaging procedure are quite different, 
while this is not for the GRBs less deviating from the red 
line. Such a naive observation makes us argue in favor of 
a problem with the measurement of the quantities entering 
the correlations of interest or of a deviation from one of the 
best fit scaling relations. Investigating in detail this issue 
on a case -by -case basis would need to retrieve the original 
data for each GRB which is outside the aim of our work. 

In order to compare the HDs from the three differ- 
ent calibration methods, we consider the values of Afi = 
l^fid{z) — iJ.{z) with fj,fid{z) the theoretically predicted dis- 
tance modulus for the fiducial ACDM model. We find : 

{Afi) = -0.15 , (A/i)„ed = 0.07 , (A/i)™» = 1.15 , 

(Afi) = -0.10 , {Afj.)med = 0.16 , (A/i)™, = 1.15 , 



(A^) = -0.24 , (A/i)™ed = -0.01 , (A^i)™, = 1.21 , 

for the F, C and LR samples, respectively. One could naively 
be surprised that the F sample (whose calibration is based 
on the same fiducial ACDM model used as reference here) 
does not give systematically smaller deviations. Actually, 
the calibration is made using only GRBs with z ^ 1.4, but 
the full sample is dominated by higher z objects so that the 
mean and median of A/i are not biased by the choice of the 
reference model. As a consequence, we find the rms Ajj, is 
roughly the same for the three samples, although the mean 
values suggest that the calibration based on the linear re- 
gression leads to fi values that are preferentially larger than 
expected. From a different point of view, higher jj, values 
at high z argue in favor of a model where the transition 
from deceleration to acceleration takes place to a larger z, 
which can be achieved (for a dark energy model with con- 
stant equation of state) pushing wo to values smaller than -1 
(i.e., going to the phantom regime) or decreasing the matter 
content. However, such a naive argument should be recon- 
sidered taking care of the uncertainties on the single GRBs 
distance modulus estimates. 

Actually, from the point of view of the impact of cal- 
ibration methods, we are not interested in minimizing the 
mean or median values of A/i, but rather on comparing the 
A/x distributions among them. Should these distributions be 
equal, then one could conclude that the HDs based on differ- 
ent calibration methods are consistent with each other. The 
different mean and median values seem to suggest that this 
is not the case. However, these distributions are quite asym- 
metric and wide enough to allow for a considerable overlap 
among them. Moreover, for each single GRB, the values of 
Ajj, are consistent with each other within the uncertainties. 
We can therefore conclude that the HDs obtained using dif- 
ferent calibration methods are consistent with each other 
within the uncertainties. 

We, nevertheless, warn the reader that such a result is 
partially due to the large error bars. Should these be re- 
duced, the consistency among the three different HDs could 
also gone lost. To this end, however, one should first reduce 
the intrinsic scatter of the correlations used in the average 
procedure since it is this term which typically dominates the 
error budget. At the present time, one should give off the 
correlations with the larger aint values (e.g, the L-V one), 
but this would lower the number of objects in the sample and 
increase the error because of averaging over a lower number 
of n estimates. Larger GRBs samples with measured redshift 
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Table 4. Deviations from the fiducial ACDM distance modulus 
estimates for the different GRBs 2D correlations used adopting 
the calibration method based on the fiducial luminosity distances 
(i.e., the F sample). We report the number of GRBs and the 
mean, median and rms values. 



Id M (Am) {/^fJ.)m.d (AA»)r. 





115 


-0.20 


0.01 


1.32 


L-Tiag 


60 


0.03 


0.09 


0.61 


L-TfiT 


79 


-0.12 


0.07 


0.94 


L-V 


115 


-0.50 


0.05 


1.86 


-^7 " ^peak 


30 


0.04 


0.28 


0.97 


^iso ~ ^peak 


100 


-0.01 


0.05 


1.01 



and observational parameters are then needed to address the 
issue of the impact of calibration methods on a safer basis. 

5.2 Impact of the averaging procedure 

As yet stated above, averaging the fi values from different 
correlations helps reducing the total uncertainties and par- 
tially washes out the systematics connected to each single 
scaling relations. However, this does not come with no price 
to pay. Indeed, such a strategy implicitly assumes that all 
the scaling relations are physical ones and not the artifact of 
a not well controlled selection effect. This danger is actually 
hard to avoid since no clearly defined and well accepted the- 
oretical model is presently available to describe the physics 
of GRBs and predict scaling relations in agreement with the 
empirically motivated ones. 

As a first check, we compare the Afj. values obtained 
estimating fi using each single correlation. For the F sample, 
we get the results summarized in TableUwhere we also give 
the number of GRBs used. The conclusions are unaltered 
for the C and LR samples so that we do not discuss them. 

While the median values of A/^ are roughly compara- 
ble, both (A/i) and {Afijrms are definitely larger for the 
L-Epsak and L-V correlations. It is noteworthy that these 
are also the only two correlations where all the GRBs can 
be used to infer the HD so that one could ascribe the larger 
(A^)rms to the inclusion of a higher number of outliers. On 
the other hand, this same effect can be a statistical artifact 
due to the use of a smaller GRBs sample. It is actually hard 
to understand whether a selection effect is at work here. In- 
deed, the inclusion of a GRB in a sample used for a given 
correlation depends mainly on observational requirements 
(e.g., the GRBs afterglow should last enough to measure 
the rise time) so that one could argue that this has nothing 
to do with its departure from a theoretical quantity such as 
the fiducial distance modulus. On the other hand, should 
selection effects induce a fake correlation, one can imagine 
that the deviations from the fiducial curve are randomly 
distributed leading to small (A^). 



Pending the questio of which relation is physical, we 
can quantify the impact of an incorrect assumption by eval- 
uating again the distance moduli excluding the L-V and L - 
Epeak correlations. Since these are the ones with the greatest 
rms AjU and the largest intrinsic scatter (see Table[l]), we ex- 
pect them to have the greatest impact on the departure from 
the fiducial AGDM theoretically predicted HD. We find : 

(A^) = -0.11 , (A/i)„ed = 0.04 , {A^l)rn.s = 1.05 , 

for 106 GRBs, being the number smaller than before because 
some GRBs with /i determined from the L-Ep^ak and/or 
L-V correlations only now drop out of the final sample. 
Compared with the case with all correlations included, the 
mean, median and rms A/^ are indeed smaller. However, 
since we average now on a lower number of correlations, the 
uncertainty on ^ for each GRB is now larger so that the A^ 
values with and without these two correlations are consistent 
within the la confidence ranges. 

We must therefore conclude that, given the available 
GRBs dataset, including or not a given correlation in the 
averaging procedure should be a compromise between the 
need to avoid an uncertain systematic bias and the desire to 
ameliorate both statistics and precision. 



5.3 Satellite dependence 

The XSIO GRBs sample is made out by collecting the data 
available in the literature so that the final catalog is not 
homogenous at all. In particular, the satellite used to get 
both the spectral and afterglows data change from one case 
to another. In order to investigate whether this could have 
any impact on the HD, we consider again the deviations from 
the fiducial ACDM model using only the 80 GRBs detected 
with the Swift satellite (excluding outliers). We get : 



(Am> = 


-0.29 , 


(A/i)med 


= 0.03 , 




= 1.33 


(Am> = 


-0.24 , 


(A/i)med 


= 0.11 , 




= 1.32 


(Am> = 


-0.40 , 


(A/i)med 


= 0.05 , 




= 1.40 



including all correlations and using the F, C and LR sam- 
ples, respectively, while we find 





-0.25 , 


(AM)med 


= -0.05 




= 1.18 


(Am> = 


-0.18 , 


(AM)med 


= 0.06 , 




1.18 , 


(Am) = 


-0.33 , 


(AM)med 


= -0.05 




= 1.24 



excluding the L - Ep^ak and L-V correlations and the same 
samples. Somewhat surprisingly, we find larger A/i values 
independent on the calibration procedure adopted and use 
of all or only four out of the six scaling laws. Actually, this 
is most likely a consequence of the different redshift range 
probed. Indeed, Afi positively correlates with z so that a 
sample with a larger median redshift will have a larger (A/i). 
This is indeed the case here since the pre - Swift data are 
limited to a smaller redshift range mainly because of instru- 
mental effects. We therefore agree with the results in XSIO 
who have argued against a systematic difference among pre 
and post - Swift GRBs datasets. 
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6 CONCLUSIONS 

GRBs have recently attracted a lot of attention as promising 
candidates to expand the Hubble diagram up to very high 
z, deep into the matter dominated era thus complementing 
SNela which are, on the contrary, excellent probes of the 
dark energy epoch. However, still much work is needed in 
order to be sure that GRBs can indeed hold this promise. 

As the Phillips law is the basic tool to standardize 
SNela, the hunt for a similar relation to be used for GRBs 
has lead to different empirically motivated 2D scaling re- 
lations. However, the lack of a local GRBs sample leads 
to the so called circularity problem, i.e. the need to know 
the cosmological model to infer the luminosity distance to 
each GRB contrasted with the desire to constrain that same 
cosmological model. In an attempt to overcome this prob- 
lem, we have here considered the impact on the scaling re- 
lations and GRBs HD of three different methods to esti- 
mate the luminosity distance. First, we rely on a fiducial 
ACDM model fitted to the SNela data. Then, we make a 
step further towards a model independent calibration us- 
ing the cosmographic expansion to the fifth order in z thus 
only assuming that the Robertson - Walker metric describe 
the universe background. Finally, we use local regression on 
SNela to interpolate the luminosity distance to a given z 
with no assumption at all on the cosmological model. We 
find that these three conceptually different methods to esti- 
mate the luminosity distance and hence calibrate the GRBs 
scaling relations lead to consistent results. Indeed, the slope, 
zeropoint and intrinsic scatter of the 2D correlations are in 
good agreement within the 68% confidence ranges. As a con- 
sequence, the Hubble diagrams averaging over the six cor- 
relations considered is not affected by the choice of the cal- 
ibration method being the distance moduli to each GRB in 
agreement within the errors. Such a preliminary conclusion 
may suggest that the calibration problem could actually be 
a false problem since the use of a fiducial model to estimate 
the luminosity distance to the z ^ 1.4 GRBs finally lead to 
the same HD of a fully model independent method such as 
the local regression. Actually, some important caveats have 
to be taken into account to not overrate this result. First, we 
have chosen as a fiducial model the concordance ACDM one 
which is known to fit very well the SNela data. As such, the 
predicted luminosity distance closely follows the one inferred 
from SNela which can be estimated using local regression. 
As such, these two apparently different luminosity distance 
methods actually give the same luminosity distance values so 
that the agreement of the constraints on the GRBs parame- 
ters is not surprising. Should we have used as fiducial model 
a different dark energy scenario, the results could have been 
different. Moreover, although not statistically meaningful, a 
weak dependence of the slope of the GRBs scaling relations 
on the calibration method is indeed present with the lo- 
cal regression leading to shallower relations. A larger GRBs 
sample is therefore needed to check whether these trends are 
only a statistically meaningless fluctuation or an evidence of 
the key role of the calibration method adopted. 

Once the calibration procedure has been adopted, one 
has still to check whether a redshift evolution of the GRBs 
scaling relations is present or not. Should indeed the scaling 
relations coefficients be a function of z, neglecting such a 
dependence would bias both the calibration and the Hubble 



diagram thus introducing a systematic error when using the 
derived HD to constrain the cosmological parameters. We 
have therefore explored two different parameterizations as- 
suming that the zeropoint only changes with z or that both 
the slope and the zeropoint may be approximated as linear 
functions of the redshift. Although there are some evidences 
in favor of a significant evolution of the zeropoint with the 
redshift, it is nevertheless worth stressing that adding one or 
two more parameters to the fit of the scaling relations does 
not improve the quality of the fit (the reduced being 
the same or larger) and do not reduce the intrinsic scatter. 
The Occam's razor therefore makes us argue that such an 
evolution is not statistically motivated and neglecting it is 
actually a conservative and better motivated choice. It is, 
however, worth saying that such a conclusion heavily relies 
on the quality of the data. Indeed, the value of the de- 
pends also on the size of the uncertainties. For instance, ar- 
tificially reducing by hand the errors leaving unchanged the 
central values increases the statistical significance of the red- 
shift evolution making this choice the most conservative one. 
On the other hand, one should reconsider this issue by pre- 
liminarily investigating whether a selection effect is present. 
Indeed, we have implicitly assumed that the probability to 
measure the quantities entering the fitted correlation is the 
same whatever the GRB redshift is. Should this not be the 
case, one has to include a prior in the likelihood function to 
model this selection effect and check whether the results are 
affected or not by the prior itself. However, such an anal- 
ysis can only be made on a case -by -case basis and asks 
for a preliminary theoretical modeling of the correlation of 
interest to simulate the full detection process. 

Assuming that no evolution is present, we have fi- 
nally checked that the derived Hubble diagrams are not af- 
fected by systematics related to the choice of the calibration 
method, the averaging procedure or the homogeneity of the 
sample. As such, the GRBs HD could be safely used as a 
tool to constrain cosmological parameters. Although these 
results are strongly encouraging, we however warn the reader 
that they are still preliminary. Indeed, what we have actu- 
ally shown is that the systematics induced by the different 
effects we have considered are smaller than the statistical 
errors due to measurement uncertainties of the GRBs ob- 
servable quantities and the intrinsic scatter of the 2D scal- 
ing relations used. It is not possible to forecast whether the 
conclusions still hold should these errors be reduced. To this 
end, increasing the number of GRBs and improving the pre- 
cision could not be the right strategy. Indeed, the error bud- 
get on the distance modulus is typically dominated by the 
intrinsic scatter so that one should rather rely on GRBs scal- 
ing relations which are as tight as possible. From this point 
of view, for instance, one should give off the L - Ep^ak and 
L-V correlations, but this comes at the price of reducing 
the number of usable GRBs and increase the error (since 
we now average on a lower number of fi estimates). Again, 
this leads to the need of larger GRBs samples to select the 
tightest correlations without degrading the precision on the 
distance modulus determination. 

As a final remark, an analogy between SNela and 
GRBs as cosmological tools can be considered. As soon 
as the Phillips law was established, the SNela started 
to be used to probe the Hubble diagram up to z ~ 1.5 
notwithstanding possible problems with the evolution of 
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the standardization method adopted and the universahty of 
its basic assumptions. Nowadays, the SNela samples have 
so many objects that it is possible not only to use them as 
cosmological tools, but also carefully explore systematics 
and their impact on the Hubble diagram. The present 
day situation for GRBs is similar with many Phillips 
law -like relations proposed to standardize them and most 
attention dedicated more to their cosmological use than 
to the systematics. As for SNela, we must therefore only 
wait that, as time goes by and new instruments bring 
us larger and higher quality samples, the issue of GRBs 
systematics can be more efficiently faced off showing that 
these powerful explosions can indeed hold the promise to 
fill the gap between the dark energy era probed by SNela 
and the early universe tested by CMBR. 

Acknowledgements VFC is supported by the Italian 
Space Agency (ASI). 



APPENDIX A 

We report here for completeness the constraints on the cal- 
ibration parameters for the two different redshift evolution 
parameterizations and the three luminosity distance esti- 
mate methods. Although the best fit values may be different 
from one case to another, the confidence ranges for the pa- 
rameters always overlap quite well so that the difference is 
not statistically meaningful. The main conclusions in Sect. 4 
qualitatively apply also changing the calibration method so 
that we do not repeat here the analysis of the results. 
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Table 1. Constraints on the calibration parameters {a,b,a,ai„t) for the 2D correlations considered in the text and fitted to the _F, C 
and LR samples (first, second and third row of each fit). Columns are as follows : 1. correlation id, 2. reduced x^i 3. best fit parameters, 
4., 5., 6., 7. median values and 68% confidence ranges for the fit parameters. The number of GRBs used in the same as in Table[T] 
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Table 2. Constraints on the calibration parameters (ao, log ai , bo > log 61, 0"int) for the 2D correlations considered in the text and fitted 
to the F, C and LR samples (first, second and third row of each fit). Columns are as follows : 1. correlation id, 2. reduced x'^i 3. best fit 
parameters, 4., 5., 6., 7., 8. median values and 68% confidence ranges for the fit parameters. Number of GRBs as in TablefTI 

Id x^/d.-0-f- (ao,logai,feo.log6i,o-inf)6/ "0 logai 60 logfei dint 
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L-V 1.24 (-0.11,-2.68,51.32,0.02,0.56) -0.2ll°;4* -2.15l37g 51.85to;53 -0.77t2;gg 0.58lo;Jg 

L-V 1.40 (-0.13,-3.13,51.40,0.06,0.43) -0.02^044 -l-94t2.39 51.86lo4g -0.5713^3 0.49115^8 

E-y-Epeak 1.70 (0.79,0.04,50.37,-0.56,0.10) lAOtoAV -^■'^^-i'sT- ^0.6lto,il -2.86j:3;g4 0.17l:o:°g 

E^-Ep^ak 1.49 (0.98,-0.05,50.46,-0.92,0.12) 1.43^0.49 -2.42+4 50.59lo;'|'2 -2.98j;3;35 0-17lo;°g 

Ej-Ep^a.k 1.70 (0.74,0.09,50.40,-0.57,0.06) 1.47l:o;3? -2.70^4.5^ 50.64l°;°* 

Eiso-Epeak 0.79 (0.94,-0.23,52.00,-0.15,0.49) 1.27^0,11 -1.781:2.54 52.56^0:34 -1.73 

Eiso-Epeak 0.91 (0.69,-0.07,51.98,-0.12,0.44) 1.20t°;g2 -l-^stsil 52.53lo;43 -1.93 



-2.41 '-'•'^^-0.08 



+1.67 0.49+"-°* 



0.07 



Eiso-Ep^ak 0.86 (0.66,-0.05,52.21,-0.28,0.45) 1.24t°-j2 -l-SSta.gg 52.57lo;3g -1.84l2:8o 0.49lo:o8 
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